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We sketch the major steps in a functional integral derivation of a new set of Stochastic Gross- 
Pitaevsky equations (GPE) for a Bose-Einstein condensate (BEC) confined to a trap at zero tem- 
perature with the averaged effects of non-condensate modes incorporated as stochastic sources. 
The closed-time-path (CTP) coarse-grained effective action (CGEA) or the equivalent influence 
functional method is particularly suitable because it can account for the full back-reaction of the 
noncondensate modes on the condensate dynamics self-consistently. The Langevin equations derived 
here containing nonlocal dissipation together with colored and multiplicative noises are useful for 
a stochastic (as distinguished from say, a kinetic) description of the nonequilibrium dynamics of a 
BEC. This short paper contains original research results not yet published anywhere. 



I. INTRODUCTION 

Two major paradigms are commonly adopted in describing the nonequilibrium dynamics of a many-body system. 
One is by way of kinetic theory based on the Boltzmann equation and its generalizations, which evolve the lowest 
order correlation functions (in the BBGKY hierarchy). The other captures the stochastic and dissipative dynamics 
of an open system, with the environment functioning as sources of noise with specific features determined from first 
principles, not put in by hand; and their overall effect engendering dissipation in the open system dynamics. Usually 
for a system well discernible from and whose effects dominate over its environment (the condensate is an open system 
with the noncondensate as its environment), the stochastic dynamics is often an easier and perhaps cleaner approach. 
Thus explains the great interest in deriving reliable stochastic equations for BEC dynamics. 

For BEC stochastic dynamics, historically Gardiner and his co-workers [3[ were the first to derive a set of stochastic 
Gross-Pitaevsky equations. Their method combines the kinetic theory and the open systems approaches. We present 
here an alternative derivation via the closed-time-path (CTP) [4J coarse-grained effective action (CGEA) [5| or the 
closely related influence functional @ method. This method has the merit that it can account for the full back- 
reaction of the noncondensate modes on the condensate dynamics self-consistently. Our results are structurally 
similar to theirs but the content is more complex: Our Langevin equations contains nonlocal dissipation and colored 
and multiplicative noises. These two kernels provide all the necessary information for a stochastic description of the 
nonequilibrium dynamics of a BEC. |8| 

II. QUANTUM FIELD THEORY OF BEC DYNAMICS 

For a field-theoretic description of BEC we begin with a second - quantized field operator ^ (x, t) which removes 
an atom at the location x at times t. It obeys canonical commutation relations 

=0, (I) 
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[*(f,t),ttt(£t)] =S(x-y). (2) 
The dynamics of this field is given by the Heisenberg equations of motion (we choose units such that U = 1) 

-i|U = [H,*], (3) 

where H is the Hamiltonian. The theory is invariant under a global phase change of the field operator 

* -» e 40 *, tft ^ e - l9 ¥. (4) 

The constant of motion associated with this invariance through Noether's theorem is the total particle number. 
We shall consider only the simplest Hamiltonian 



H = J d d x j*t#tf + ^1^2 j 



(5) 



where we have reduced the atom - atom interactions to a single contact potential with a delta-like repulsion term. 
The coupling constant U is often parameterized in terms of the scattering length a and the mass M of the atom as 
U = Aita/M. The single - particle Hamiltonian H is given by 

Hy = -^VH + V trap (x)V, (6) 

Vtrap (x) denotes a confining trap potential. Then the Heisenberg equation of motion 

d 

dt ' w 

is also the classical equation of motion derived from the action 

S = J dt d d x i**-^* - J dt H. (8) 



III. CONDENSATE AS OPEN SYSTEM 



We adopt the quantum open system framework and define our system of condensate as comprising a superposition 
of only a few low lying modes of the one-particle Hamiltonian. The condensate while evolving can excite the quantum 
fluctuations in the higher modes, which in turn back-react on the condensate, thus modifying its test-field dynamics 
(i.e, that which is obtained without backreaction) . We regard the condensate as an open system living in the 
environment provided by the higher modes and try to incorporate the effects of the coarse-grained environment in 
some way. The inclusion of backreaction effects in our experience is best performed via the coarse grained effective 
action or the influence functional as they ensure full self-consistency. 

We first partition the total wave function into the system sector and an environment sector. We do this by way of 
projection operators. Assume a complete set of one-particle states 

Htp a (x) = uj a tp a (x) . (9) 

Choose an appropriate partition frequency uoq and define two projectors 

p= iM2)iM20> ( 10 ) 

which projects onto the low lying modes (band) comprising the condensate, our system, and 

Q = 5(x-x')-V. (11) 

which projects onto the sector of higher modes, the noncondensate, constituting an environment for the system. 
Define the condensate band field operator 



4> (t, x) = TV (t, x) , 



(12) 
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and split 

= 4>(t,x)+ X (t,S). (13) 

Here, (j) denotes the quantum wave function of the condensate, not a c-number wave function usually assumed. 
We now split the action into condensate (C), non-condensate (NC) and interaction (int) parts 

S [<£ + X] - S c [0] + S NC [x] + S int [0, X] • (14) 

The first two are just the full action evaluated at the corresponding field. Sj nt may be written as the sum of three 
terms in successive higher powers of the noncondensate operator: 

Sint = Si + S 2 + S 3 , (15) 

with 

Si = -U J dt d d x {^fa + xW} , (16) 

S 2 = -| J dt d d X {^V + x t2 2 + 400x t x } ; (17) 



S 3 = -U J dt d d X {0t x t x 2 + x t2 x0 } ; 



(18) 



where eft is the hermitian conjugated field of <j>. To simplify the appearances we represent the pair (0, </>t) as the up 
and down components (respectively) of a spinor a , a = 1,2, and adopt the DcWitt convention where the space-time 
arguments are assumed to be included in the spinor indices. 

IV. CGEA AND STOCHASTIC GROSS-PITAEVSKY (LANGEVIN) EQUATION 

We now introduce the coarse-grained effective action. Assume for simplicity that at t = the quantum state of the 
gas is described by a factorizable density matrix 

P [4> 1 + x\4> 2 + x 2 A =pc[^ 1 ,^ 2 ,o]p JV c[x 1 ,x 2 ,o]. (19) 

The expectation values of field operators may be derived from the usual closed time-path (CTP) generating functional 
exp{iW[j\j 2 ]} = /w'^eW-^'^/^'^'-^DptO). (20) 

However, if we are interested in the expectation values of the condensate band operators, then we only need to couple 
sources to these fields, namely 

QJ 1 ' 2 = 0. (21) 
The integral over non-condensate fields may be performed, and we obtain the condensate generating functional 

exp{iW c [J\J 2 ]} = j DcfDcf e'^^^+J^' [W-^]) pc (Q) ^ (22) 

where the coarse-grained effective action is 

Scgea [<?, <?] = S c - S c [<f> 2 ] + S IF [<f>\ ^ 2 ] , (23) 
where Sif is the influence action and 

exp {iS IF [ct>\4> 2 } } = j Dx 1 D X 2 e i (^o[x 1 ]+s int [0 1 ,x 1 ]-S NC [ x 2 ]-S int [^,x 2 ]) /9jvc (0 ) . ( 24 ) 



4 



is the influence functional (IF). 

To make the physical content of the CGEA more explicit, it is convenient to introduce new field variables 

^+ = ^(^+0 2 ), <t>- = <t> 1 -<f> 2 . (25) 

If we Taylor expand Scgea hi powers of 4>- we find that there is no independent term, and that even terms are 
imaginary and odd terms are real. The path integral admits a saddle point with respect to 0_ at 0_ = 0. We 
concentrate on the contribution of field configurations near this saddle point by keeping only quadratic terms in 
in the CGEA 

S IF [<j>+, 4>-\ = J dtd d xV a [4>+]^_ 

+ l -Jdtd d xJ dt' d d x" Nab [0+] 0- 4>- + ■ ■ ■ (26) 
To get an explicit representation, we introduce the functional expectation value 

{A} [(/>+} = J Dx 1 D X 2 e i ( Sjvc [ xl ] +5 -*K^ 1 ]- s ~-o[x 2 ]-Sin*[0+,x 2 ])^ pjvc ( ) , ( 27 ) 

Then 

v - ™ - If + 5 < W [0t ' x,] + W [0+ ' x2] ) 10+1 ' (28) 

From here one can deduce that generally a representation 

M ab [0+] = ]T ^;,,,,/'; [<!>+] ft [<!>+] > (30) 

can be obtained in terms of a set of local functionals /,f . This allows us to rewrite the condensate generating functional 
as 

exp{iW c [J + ,J-}} = J DC ac T[C ac ] Dj+D^e'S 



and 



where the Feynman- Vernon trick of using a Gaussian functional identity [6J involving the stochastic distribution T is 
invoked. Here 

J+ = l(J 1 + J 2 ), J- = J'-J 2 , (32) 
and the are the Gaussian stochastic sources with zero mean and cross correlation 

£ac£bd) ~ u ac,bd- (33) 



It should be noted that the stochastic terms do not appear explicitly, but only implicitly in the form of random 
contribution to the action, which then is averaged over. However, performing the final integration over § a _ we see 
that the only contributions to the path integral come from condensate field configurations which obey the equation 



v\v a [0+] + 2 e ac f c i [</>+] + j+j = o, 



(34) 
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with random initial conditions weighted by the Wigner function constructed from the initial density matrix pc (0). 
The free evolution of the condensate is obtained by setting J + = 0. In this sense, the condensate field evolves as 
an ensemble of trajectories, each obeying the stochastic Langevin-type equation (|34|) . This is a consequence of the 
truncation of Scgea to just quadratic terms in <f> a _. Further terms in Scgea would introduce effects associated to the 
quantum nature of the condensate field. The presence of the projection operator V in (|34[) enforces the consistency 
of the system - environment split under time evolution. 

We observe that because of the nonlocal terms in D a this equation is generally dissipative and contain noises 
which are generally colored because of the nonlocality of the cross correlations v? bd , and multiplicative because of 
nonlinearity in the ff. There ought to be a fluctuation-dissipation relation for each types of the noise. A similar 
description for a theory with polynomial type of system-environment coupling can be found in Hu, Paz and Zhang 
(1993) of @|; see also Q. 

V. DISSIPATION AND NOISE KERNELS 



We can glean more physical meanings about dissipation and fluctuations by trying to obtain more explicit expres- 
sions for these kernels. We go back to our model of a gas confined in a time-independent trap, and compute the 
CGEA to order U 2 at zero temperature (T = 0) for the vacuum. We obtain (from here on, to simplify appearances, 
we drop the + subindex on <f> + ) 

a M = + — J dt dx 

+ i^(t,f)^ +2 (t-i',x,x')0 2 (t',f) 

~ 1%:^)^^-^^)^^',^, (35) 

where 

F+(t-t>,x,x>) = (x(t,x) X Ht\x')) 

= e-^'-'V^)^'), (36) 

Wot >^C 

F- (t-t',x,x') = ( X {t\x')x ] {t^)) 

= e^-^^W^iaf). (37) 

u a >ui c 

With respect to the noise terms, we find that they may be represented in terms of four stochastic terms C 1 to £ 4 

& e = 8 ac C i , i = l,... ,4, (38) 



coupled to four functions /i to /, 
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h« = 4£;W 2 <f>-<t> i <t> 2 }, (40) 



2 d<j) a 



/- = ^[^+^]. ( 41 ) 



The nontrivial cross correlations are 



where 



(C 1 (t, x) C 1 (/, S> )) = <C 2 (t, x) C 2 (f , £)) = \f p (t - t', x, f ) , (43) 

(C 1 (t, S) C 2 (*', f )) = ^ (i - f, f, f ) , (44) 

<C 3 (t, a?) C 3 (f , f )) = <C 4 (*, 5) C 4 it', x' )) =F p (t- f, x, x") , (45) 

<C 3 (t, x) C 4 (f, x')) = iF m (t - t', x, x') , (46) 

F p =F++F-, F m = F + — F~, F p = F+ 2 +F~ 2 , F m = F +2 - F~ 2 . (47) 



The detailed forms of the noise and dissipation kernels in the stochastic Gross-Pitaevsky (Langevin) equation are 
useful for a comprehensive and in-depth study of the stochastic dynamics of the condensate incorporating its inter- 
action with the noncondensates. Various aspects of this problem are understudy and results will be reported elsewhere. 
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